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WSTĘP 


W poprzedniej pracy [1] przedstawiono metodę symulacji cyfrowej do ba- 
dania procesu EC I (reakcja chemiczna pierwszego rzędu poprzedzona procesem 
elektrodowym) na wirującej elektrodzie dysk-pierścień oraz omówiono zasadnicze 
części napisanego w Fortranie programu do symulacji tego procesu. Obecnie zaj- 
miemy się rozszerzeniem powyższego programu tak, aby był on przydatny do ba- 
dania kinetyki reakcji homogenicznej drugiego rsędu w warunkach miareczkowania 
w warstwie dyfuzyjnej (proces EC II). 

Od strony analitycznej problem opisu procesu EC II sostał w sposób przy- 
bliśony roswiąsany przez Albery'ego i in. [2, 3, AL W pracy [4] podano sposób 
wyznaczania stałej szybkości ką reakcji drugiego rsędu z krzywych zależności tę 
od ip dla różnych szybkości wirowania elektrody oras dla różnych stężeń substan- 
cji C reagującej s substancją B generowaną na dysku. W skrócie metoda polega 
na wysnaczeniu tsw. kinetycznego prądu dysku px, tj. wartości prądu dysku w 
momencie gdy powierschnia reakcji substancji B i C snajduje się na wewnętrznej 
krawędzi pierścienia [3,4] (powierzchnią reakcji określa się obszar, gdzie [B] = ICH. 
Przy takim prądzie dysku przy stałej prędkości wirowania elektrody w (dla da- 
nego stężenia substancji C) prąd pierścienia tp jest uzależniony od ilości substancji 
B przenikającej w obszar dominacji substancji C, a więc w tych warunkach ip 
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zależy od szybkości reakcji pomiędzy B i C. Mając wyznaczone wartości (Gr dla 
różnych w i dla różnych stężeń substancji C, określa się wartości kinetycznego prądu 
pierścienia (Sr s krzywych tp vs. iR. Mając te wartości można określić zależność 
Ny od w/C$, gdzie Ny = |ip.k/tp,k| a CQ jest stężeniem substancji C. Zgodnie z 
teorią Albery'ego i in. [4] Nu wyraża się wzorem: 


Np = .339r3r; DV 3y"1/8|1 — F(a)|wkaC$ (1) 


gdzie r; i r2 oznaczają promień dysku i promień dysku s przerwą izolacyjną, D 
jest współczynnikiem dyfuzji, v — lepkością kinematyczną rostworu a F(a) jest 
funkcją wymiarów elektrody [5]. Widać, że wykres zależności Nr od w/CQ, powinien 
być linią prostą przechodzącą przez początek układu. Jest to prawdą jedynie dla 
bardzo małych wartości w/C%. Dla większych wartości w/C wartość N; nie rośnie 
nieograniczenie (jak to wynika ze wzoru) lecz asymptotycznie dąży do wartości No 
[4,5] (efektywności zbierania elektrody). Wyznaczenie ką z zależności Nr vs. w/CQ 
polega na określeniu nachylenia prostoliniowego odcinka tej zależności, które jest 
proporcjonalne do 1/ką. , 

Prater i Bard przedstawili metodę wyznaczania stałej szybkości ką opartą 
na symulacji cyfrowej [6]. Metoda polega na symulacyjnym wyznaczeniu krzywych 
roboczych Nz od vim (gdzie x? = kąC93w”1u1/3D-1/3(0.51)-2/3, m = C3,/C9, 
CY oznacza początkowe stężenie substancji elektroaktywnej) przy skoku potencjału 
na dysku do wartości prądu granicznego. Otrzymuje się w ten sposób (dla danej 
elektrody) rodzinę krzywych — dla każdej wartości parametru m jedna krzywa. 
Ustalając w eksperymencie wartość «p na wartość graniczną tp, (graniczny prąd 
dyfuzyjno-konwekcyjny) należy rejestrować zmiany NM (= |tRu/tp„|) przy smianie 
w i CQ. Wartość ką można otrzymać przes wyznaczenie s krzywych roboczych 
wartości x?m dla każdego Nu i wykreśleniu zależności vim od w/C%. Nachylenie 
prostej jest wprost proporcjonalne do 1/k2. 

Naszym celem jest zaprezentowanie innego sposobu wykorzystania symulacji 
cyfrowej do wyznaczania ką. W symulacji przez nas proponowanej na dysk nakłada 
się prąd o wartości takiej, aby powierzchnia reakcji znajdowała się na wewnętrznej 
krawędzi pierścienia. Dla otrzymanych z analizy doświadczalnych krzywych ip 
vs. tp wartości Nr |4] znajdujemy odpowiadające im wartości vim na roboczej 
krzywej symulacyjnej Ną od vim. Następnie wykreślamy zależność vim od Cou 
(zależność prostoliniowa). Nachylenie prostej jest wprost proporcjonalne do 1/ką. 

Metoda ta ma kilka zalet w stosunku do metody Albery'ego i in. Między 
innymi charakteryzuje się większą dokładnością i jest użyteczna bes nakładania 
ograniczeń na wartości w/CQ. Zaletą jej jest także to, że dla danej elektrody 
mamy jedną krzywą roboczą Nr vs. vim bes wsględu na wartość parametru m. 
U Pratera i Barda dla każdego parametru m należy wykonać osobne obliczenia. 
Prócz tego, jak wskazuje Albery i in. [4] dla wyznaczania ką najkorzystniej jest 
„ustalić” pozycję powierzchni reakcji na wewnętrznej krawędzi pierścienia, gdyż 
wtedy pierścień „sbiera” jedynie tę część substancji B, która nie sdążyła przerea- 
gować s substancją C w rejonie jej dominacji. Tak więc opisywana metoda łączy se 
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sobą zaletę metody Albery'ego i in. (praca w zakresie prądów kinetycznych) oras 
zaletę metod symulacyjnych (brak ograniczeń sakresu stosowalności). 


SYMULACJA CYFROWA PROCESU EC II ZE SKOKIEM PRĄDU NA 
DYSKU 


Symulowany proces można opisać równaniami: 


A+tne  B dysk (2) 
B+C€—D+E warstwa dyfuzyjna (3) 
B+tne A pierścień (4) 


gdzie zakłada się, że substancje B, C i D są elektronieaktywne. W stosunku do 
procesu EC I powyższy proces różni się jedynie reakcją (3). Do symulacji opisa- 
nej poprzednio należy kilka zmian związanych z II-rsędowością reakcji w warstwie 
dyfuzyjnej, skokiem prądu na dysku oraz dyfuzją i konwekcją substancji C. Wpro- 
wadsamy parametr m, określony zależnością [6]: 


nm C/O 


gdzie C9, i C9, oznaczają początkowe stężenia substancji C i A (stężenia substancji 
A i C w głębi rostworu). Początkowe stężenie substancji C w każdej komórce siatki 
przestrzennej będzie równe [1]: 


FC1(J, K) =m 


Dyfuzję i konwekcję substancji C uwzględniamy analogicznie jak dyfuzję i 
konwekcję substancji A i B [1, 6]. Zmianie musi ulec parametr zwiąsany s szybkością 
reakcji homogenicznej. Dla reakcji II rzędu mamy: 


—dCo/dt = —dCp /dt = kaCzCo 


Wyrażając to przy pomocy parametrów symulacyjnych (tsn. Cc/C3 = FC, 
CB/C3 = FB) otrzymamy [1, 6]: 


—AFC = -AFB = ką FB CcAt 
Ponieważ D, 7] At = E = v!/8D-1/5,-1[-1(0.51)-2/8 a Co = FC CY, 
to -AFC = -AFB = kąt4C3FB FC/L. Beswymiarowy parametr Eat? 
ozsnaczmy XKTC [6]. Mamy: 
XKTC = bat? = kąC9w”1v!/3DZ*/*(0.51)72/3 


A więc ostatecznie 


FB1(J, K) = FBi(J, K)(1— XKTC FC1(J,K)/L) 
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FC1(J, K) = FC1(J, K)(1— XKTC FB1(J, K)/L) 


Ostatnia smiana dotycsy warunków brzegowych przy skoku prądu na dysku. Na 
dysku mamy: 


FA1(1, 1) = FA1(1,1) — FLUX — DMA(FAi(1, 1) — FA1(2,1)) 
FB1(1,1) = FB1(1,1) + FLUX — DMB(FBi(1, 1) — F B1(2, 1)) 


FC1(1, 1) = FC1(1, 1) - DMC(FC1(1, 1) — FC1(2, 1)) 


W przerwie isolacyjnej i na pierścieniu warunki brzegowe są takie jak po- 
przednio [1] (oczywiście należy uwzględnić dyfuzję substancji C do elektrody). 
Parametr F LU X określa ilość substancji A jaka uległa reakcji elektrodowej 
na dysku w czasie At przy założeniu, że na dysk nałożony jest prąd „kinetyczny” 
Ponieważ ogólnie [7] 


FLUX = DMA(FA1(2, 1) — FA1(1, 1)) = ZD(DMA/L)'/2 
a [1, 7] 
ZD = ip /(0.51)!/ nF AC9 DĄw!/2y-1/8 


(ZD jest prądem dysku wyrażonym przes wielkości beswymiarowe) 
oras . 


1 = .62n FACH D2/3w1/2y-1/6 


(*. jest prądem dysku w stanie stacjonarnym) 
to możemy napisać 


ADp/i, = (0.51)*/*/0.62FLU X(L/DMA)'/2 


czyli 
FLUX = (ip/i,,)0.62/(0.51)1/*(DMA/L)V/2 


Ponieważ śp musi być równy prądowi „kinetycznemu” dysku, który spełnia 
zależność [2, 5]: 


Kp,x| = M/(1 — F(a)) 


gdzie M = nFAD2/*(0.51)1/3w1/2y-1/6C9, /1.288 a F(a) jest funkcją saleśną od 
wymiarów elektrody, przy czym a = (ra/r1)3 — 1, to ostatecznie możemy sapisać : 


FLUX = |m/(1 — F(a)]0.776(DMA/L)/2 
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XKTC:-m 


Ryc. 1. Wykres zależności Ne od XKTCm dla elektrody o wymiarach rı = 0.3 cm, 
= 0.33 cm, ra = 0.42 cm. o m = 0.00871, [ m = 0.004354, A m =0.006531, x 
m = 0.001088. 


WYNIKI I DYSKUSJA 


Dla danej elektrody wartość Nz jest taka sama dla układów s tą samą 
wartością vim, przy csym prąd dysku musi być tak dobrany, aby powierzchnia re- 
akcji snajdowała się w tej samej pozycji dla kaśdego s tych układów. Innymi słowy 
układy s tą samą wartością x?m i tym samym położeniem powierschni reakcji dają 
te same wartości Nr. Dsięki temu w symulacji procesu EC II prsy skoku prądu 
na dysku do wartości, prsy której powierschnia reakcji osiąga wewnętrsną krawędź 
pierścienia, otrzymujemy jedną krsywą roboczą saleśności N, od XKTC m, bes 
wsględu na wartość m. Krsywą taką przedstawia ryc.1 . 

Opisaną tutaj metodę użyto do wysnaczenia stałej szybkości reakcji roda- 
nowania fioroglucyny. Na ryc.2 przedstawiono zależność X KTC m od C%/w, przy 
czym wartości X KTC m odpowiadają doświadczalnie wysnacsonym wartościom 
Ny. 

Rycina 3 przedstawia eksperymentalną krsywą N od w/CQ. Zgodnie s 
sależinością Albery'ego, początkowa prostoliniowa część tej krsywej odpowiada 
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Ryc. 2. Wykres saleśności XKTC m od C$/w [mol s/dm*rad] g C8 = 1.3 107% 
[mol/dm*], (0) C$ = 1.7 107* [mol/dm*|, (A) C8 = 8.7 1074 [mol/dm*] 


saleśności (1). W pracy [8] porównano teorię Albery'ego i in. s metodą symula- 
cji cyfrowej, sugerując smianę współczynnika .339 we wsorse (1) na .424, dsięki 
czemu wyniki symulacji są bardsiej sgodne s teorią. Analisując ryc. 2 i 3 widać 
więkssy sakres oras większą precysję metody symulacji cyfrowej . 
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Ryc. 3. Wykres zależności Nr od w/C& [rad dm*/mol a), -] C3 = 1.3 107° |mol/dm*|, 
(A) C8 = 1.7 107°? [mol/dm?], (o) C8 = 8.7107* |mol/dm*] 


LITERATURA 


DI Penar J., Sykut K., Annales Unie. Mariae Curie-Sklodowska, 43, 44, (1988/1989) 

[2] Albery W.J., Bruckenstein S., Johnson D.C.; Trans. Faraday Soc., 62, 1938, (1966). 
Il Albery W.J., Bruckenstein S.; ibid., 62, 2584, (1966). 

[4] Albery W.J., Hitchman M.L., Ulstrup .; sbid., 65, 1101, (1969). 

[5] Albery W.J., Hitchman M.L.; Ring-disc electrodes, Claredon Press, Oxford 1971. 
[6] Prater K.B., Bard A.J.; J. Electrochem. Soc., 117, (3), 335, (1970). 

[7] Prater K.B., Bard A.J.; ibid.,117, (2), 207, (1970). 

[8] Albery W.J., Drury J.S; J.Chem.Soc.Faraday Trans., 68, (3), (1972). 


SUMMARY 


This work presents a new way of usage the digital simulation for second order kinetics 


study on rotating ring-disc electrode. 
Złożone w Redakcji 23 II 1990 
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